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When considering regularity of surfaces, it is its geometry that is of interest. Thus, the con- 
cept of geometric regularity or geometric continuity of a specific order is a relevant concept. In 
this paper we discuss necessary and sufficient conditions for a 4-patch surface to be geometrically 
■ continuous of order one and two or, in other words, being tangent plane continuous and curvature 

continuous respectively. The focus is on the regularity at the point where the four patches meet 
and the compatibility conditions that must appear in this case. In this article the compatibility 
conditions are proved to be independent of the patch parametrization, i.e., the compatibility con- 
ditions are universal. In the end of the paper these results are applied to a specific parametrization 
such as Bezier representation in order to generalize a 4-patch surface result by Sarraga. 
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In many applications in Computer-Aided Geometric Design (CAGD) and Computer Graphics a sur- 
face is composed of several patches, where a patch usually is represented by a Bezier polynomial, 
B-spline or NURBS. In particular, each patch is as regular as is needed. Thus, when considering 
smoothness of a surface such as tangent plane continuity or curvature continuity, the lack of regularity 
may only occur somewhere at a common boundary curve between two or more patches. 

Regularity for a surface constituting of two adjacent patches sharing the same common boundary 
curve, see Figure 1, is a well studied problem. General results for such a 2-patch surface in the case of 
G 1 as well as G 2 continuity, i.e., tangent plane continuity and curvature continuity respectively, were 
given by Juergen Kahmann in a paper from 1983, see Q. In the same paper he applies these results to 
the case of Bezier patches. Other authors such as Degen [3], Liu and Hoschek ifTUl . Liu (9l, DeRose 
[4 ], have also treated tangent plane continuity in the 2-patch case. In the case of curvature continuity 
of 2-patch surface we refer to articles by Kiciak (H, Ye, Liang and Nowacki lf20l . 

A more complicated situation is regularity of a surface consisting of four patches where every 
pair of adjacent patches meet at a common boundary curve and all the patches intersect at a common 
vertex, see Figure |2] Among the many authors that have treated regularity problem in this 4-patch 
surface case are Bezier [ Q, Sarraga fl6l . ifTTl and [ 18 ], Ye and Nowacki fl2~T1 . Further references and 
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an overview can be found in the book by Hoschek and Lasser [6]. In refered articles so far there have 
only been considered certain cases of parametric patches, not a general parametrization as is done in 
this paper. A general approach has also been done by Peters, see the articles by Peters [L2l- lfl5l and 
by Ye Q5J. 

A general approach to study regularity for a 4-patch surface has been to restrict the patch parametri- 
zation to a certain explicit polynomial or rational basis. In this paper we consider 4-patch surfaces 
where the patches are given by any function of the form (u,v) \-¥ r(u,v) £ I 3 , u,v € [0, 1]. We 
present necessary and sufficient compatibility conditions in order to have tangent plane continuity and 
curvature continuity respectively for 'such surfaces. The results we achieve in this paper are indepen- 
dent of the patch parametrization. Thus, these results are of a general nature and can be applied to any 
parametrization. 

In the last section of this paper we apply our results by considering the Sarraga case of filling a 
hole of a G 1 -surface in such a way that the extended surface preserves the G l -regularity, see Sarraga 
lfl6l . Compared to Sarragas result we reduce the bi-degree to (5,5) of the created interior patch. More 
generally, we give necessary and sufficient conditions in order to create such an interior patch. 

2 Geometric continuity of order 1 

When discussing regularity of a surface our focus is on the geometry of the surface and not on its 
actual parametrization. Thus the notation of geometric continuity is the concept used in this context. 
The lowest order of regularity is G°, which means that the surface is connected. Another way to 
put it is to require that its representation is continuous. The next level of regularity is tangent plane 
continuity, denoted by G 1 , which is defined here. 

Definition 1 A continuous surface is said to be tangent plane continuous, denoted by G , if every 
point on the surface has a unique tangent plane, which varies continuously on the surface. Such a 
surface is also said to be geometrically continuous of order one. 

Consider Definition Q] in the case of a 2-patch surface S. Here we use the notation r € for 
a patch described by a continuous differential function (u, v) h-> r(u, v) G M 3 with < u, v < 1 
satisfying r u x r„ ^ for < u, v < 1. Let the two patches be described as (u, v) (-)• rW (u, v) with 
< u, v < 1, and (s,t) i-> A 2 \s,t) with < s,t < 1. Suppose further that each patch is regular 
enough, i.e., r^,!^ 2 ) € C^. In order for the surface S to be tangent plane continuous the only 
points that do not automatically fulfill the G 1 -condition are those along the common boundary curve 
of the two patches, see Figure 1. On this boundary curve we must particularly have v i-> rW(l, v) = 
r-( 2 )(0, t{v)) for < v < 1, where v h-> t(v) is a regular reparametrization of the interval [0, 1]. 

Thus, at a particular boundary point with parameter value v the tangent plane must fulfill 

span{r£\l, v), r« (1, v)} = span{r^ (0, t(v)), vf } (0, *(«))}. 

Let us formulate this in an alternative way, where we use the notation r-( 2 )(ti,w) = r^ 2 ) (u, t(v)) for 
< u, v < 1. The statement can be found in e. g. [20] and is summarized in the next Lemma. 

Lemma 1 A necessary and sufficient condition for two adjacent C^-patches and joining 

G 1 -continuously along its common boundary curve v i— > r^'(l,v) = r^ 2 \0,v) is that there exist 
continuous functions Ai_2 and K\% such that 

r£\0,v) = \ li2 (v)r^(l,v)+K h2 (v)r^(l 1 v), < v < 1. (2.1) 



Figure 1 : Two patches connected by a common boundary curve 



In this paper we will consider the problem of a surface constituting of four patches where every 
two adjacent patches have a common boundary curve, see Figure |2] Moreover, the four patches 
intersect at a common vertex V. In the case of a 4-patch surface being G 1 we will prove that there 
must exist compatibility conditions at the intersection point V. The compatibility conditions can be 
rephrased in such a way that we formulate necessary and sufficient conditions on the functions Xij 
and Kij at the intersection point V. 

First, by using the relation (12.11 ) we get the next four relations between the patches (l)-(2), (2)- 
(3), (4)-(3) and (l)-(4). We use the same parameters u and v for all the patches, where < u, v < 1. 
Thus 



X 1>2 (v)r^(l,v)+K 1>2 (v)r^(l,v) 
X 4)3 (v)r^(l,v) + K 4>3 (v)rW(l,v) 



(2.2) 



and 




(2.3) 




Figure 2: Four patches connected in a common vertex V 



Considering a tangent plane continuous 4-patch surface, it must satisfy that every pair of its 
patches with a common boundary curve coincide along that curve. The same must be true of its 
derivatives along the same boundary curve. 
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Thus 



(2.4) 

ri i] (l,v) =r ( v\0,v), v G [0,1] 

and 

ri 1) (n,l)=ri 4) ( U ,0) 

ri 2) (n,l) =ri 3) (u,0), nG [0,1]. 

Using the equations (12- 2b — (12.5b . we prove the next theorem, which is the main result in this section. 
In the case where the patches are described by polynomials, this result has already been published by 
Bezier in 1986. See Iffl. p 44-46. Here we formulate and prove the result for a general parametrization. 

Theorem 1 Let S be a 4-patch surface consisting of C^-patches i = 1, . . . , 4. Let V be the 
intersection point of the four patches. Then, necessary and sufficient conditions in order for the 
surface S to be tangent plane continuous are that there exist continuous functions Xij and Kij, i = 
l,j = 2,4 and i = 2,4,j = 3, satisfying the equations \2.2\ and ri2.JP , and the following relations 

«l,2 = ^1,4 ^4,3 

(2.6) 

«1,4 = M,2 «2,3 

and 

Al,2 - A 4 ,3 = «1,4 «4,3 ^ 7 ^ 

Al,4 — A2,3 = «1,2 ^2,3 

at the vertex V. 

Remark. The notations K\p, Ai ; 4, ^3, etc, are to be interpreted as ^2(1), Ai j 4(l), ^4,3(0), etc. i.e., 
as the value of the functions at the vertex V. 



Proof. In order to prove the above statement, we must see under what conditions the equations (12.21 ) 
and (12.31 ) are all satisfied at the vertex V. We use here the short notation r^ 1 ^ for r„ (u,v) |y = 
rW(l,l),etc. 

We start by eliminating ri , , r^ 4 ^ and r« in the equations (12.21 ) and (12.31 ) by using (12.41 ) and 

-ri 3 > =0 



231 ) to get 



M,2*u 


, K _(1) 


M,3 r u 




(1) 






*2,3 r v 



r u + K4 ) 3r„ — U 
-r® =0 
^ -r® =0. 



(2.8) 



With the use of (12.81 ) we replace the second and fourth equation in the above system by eliminating 
rl 3 ^ and to get 

(A 4 ,3 - Al, 2 + «l,4«4,3)r« + (Al,4«4,3 _ «i,2)r« = 
(Al,2«2,3 - Kl^TiP + (A 2 ,3 - Ai )4 + Kl ) 2«2,3)r«^ = 0. 
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Since the vectors r„ , r{, span the tangent plane, it follows from the above equations that (12.61 ) and 
(12.71 ) must hold. This concludes the proof. □ 

Let us look at some simple consequences of Theorem 1. Obviously, the functions Ajj are not 
allowed to be zero if tangential continuity is to be satisfied. Thus, if e.g. ^4^(0) ^ then /ci,2(l) ^ 0, 
which follows from (12.6I ). On the other hand, if ki,4(1) = then also ^2,3(0) = 0. This situation 
is exemplified in Figure 3. In general, it follows from equation (12.61 ) that the pair of rey's in each 
equality must both be zero or non-zero. 




Figure 3: Four patches connected in a common vertex V with «i,4(l) = ^2,3(0) = and 

Kl,2(l)«4^(0) + 

Another observation that can be done from (12.61 ) and (12.71 ) is that the next relations are true 

^1,4^4,3 = Ai j 4(Ai ! 2 — «1,4«4,3) = ^1,2^1,4 ~ Ai^K^Ki^ 

= Ai^Ai^ — ki,2«i,4 

and 

Al,2A2,3 = Ai j 2(Ai j 4 — ^1,2^2,3) = Ai^Ai^ — Ki j 2Ai i 2^2,3 

= Ai^Ai^ — ^1,2^1,4. 

In particular, we have 

Al,2A2,3 = Ai j 4A4 j 3, (2.9) 

which will be useful in the next section. 

3 Geometric continuity of order 2 

In this section we consider 4-patch surfaces of higher regularity. Therefore we introduce the following 
concept. 

Definition 2 A tangent plane continuous, G , surface is said to be curvature continuous, denoted by 
G 2 , if every point on the surface has a unique Dupin indicatrix, which varies continuously on the 
surface. Such a surface is also said to be geometrically continuous of order two. 
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Another equivalent way to describe the notation of curvature continuity is to say that the normal 
curvature at each point and in each tangential direction^ has to be unique, or the principal curvatures 
are unique. These differential geometric notations are introduced and explained in any book about 
differential geometry, e.g. GJ. First we consider a 2-patch surface, where the patches are of regularity 
C\ = C\ n C 2 . Thus, lack of G 2 -regularity for the surface S can only occur at the common 
boundary of the two patches. A necessary and sufficient condition for a 2-patch surface to be curvature 
continuous can be found in a paper by Juergen Kahmann [7]. A proof of this well-known result can 
also be found in GUI. The result is summarized in the next Lemma. 



Lemma 2 Let S be a 2-patch surface consisting of the C^-patches (u, v) rW (u, v) and (u, v) h > 
r^ 2 \u,v), < u, v < 1, satisfying equation d2.il ) along their common boundary curve. A necessary 
and sufficient condition for the surface S to be curvature continuous is that the following relation is 
fulfilled 

r£)( ,t,) = A?»r$(M +2\ h2 (v) Klj2 (v)rW(l,v) + «? i2 («)r«(l,t;) 
+ p h2 {v)r^(l,v) + iy 1 ^(v)r ( p(l,v), < v < 1, 

where the functions Ai )2 , Ki, 2 , 1*1,2 and v\fi are continuous. 

We now continue to consider a 4-patch surface as shown in Figure [2] The first and most obvious 
conditions to have geometric continuity of order two are that the equations in (12.21 ) and (12.31 ) are still 
satisfied when differentiated with respect to the parameters u and v respectively. We assume here that 
the functions Ay and Kij are differentiable. Thus 



r<$(p,v) = X' h2 (v)rW(l,v) + \ 1>2 (v)t£>(1,v) + k' 1>2 (v)v^ (1,v) + K h2 (v)r$ (1,v) 
and 

r( 4 n } (n,0) = Ai i4 («)rW(t*,l) + Ai )4 («)r«(«, 1) + k' m (u) r«(«, 1) + K M («)r«(«,l) 
r$(u,0) = A' 3 ( U )r( 2 )(n,l) + A 2 , 3 («) r^(«, 1) + k'(u)v^ (u,l) + n 2 , 3 (u)r^( U ,l) 



(3.11) 



(3.12) 



must hold. From Lemma |2l we know that curvature continuity implies that the next four relations 
must be fulfilled 

r«(0,«) = A 2 »r«(M) + 2A li2 ( U ) Kli2 ( V )r( 1 J(l, ? ;) + K?>)rW(M) ^ 
+ /i 1)2 (u)rW (1, u) + i/i, 2 (v)r( 1 )(l, v) 

r$(0,«) = Xl 3 (v)r^(l,v) + 2X i , 3 (v) K ^ 3 (v)T^(l,v) + ^(v)!® (l,v) 



and 



r®(«,0) = Al )3 (u)r$(u, 1) + 2\ 2;3 (u)k 2 ^u)t^(u, 1) + «i )3 («)^(«, 1) t 
+ ^2,3 (-")r(, 2 )(u, 1) + i/ 2 ,3 (n)rl 2 )(n, 1) 



In fact, it is enough that it holds for 3 pairwise linearly independent tangential directions by the 3-Tangent Theorem, 
see Pegna & Wolter 1111 or Hoschek & Lasser |6|, p 333. 
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r<£(u,0) = A? i4 («)r«(«, 1) + 2A M (u)K M (u)rM(u, 1) + K 2 , 4 (u)r«(u, 1) 

+ /i ll4 (u)rW(u J l) + »i ) 4(«)rW(u } l). 

The main result of this section is a general result. It concerns parameterized patches of any kind and 
it is formulated next. 

Theorem 2 Let S be a G 1 4-patch surface consisting of C^-patches r™, i = 1, . . . , 4. Let V be 
the intersection point of the four patches. Then, necessary and sufficient conditions in order for the 
surface S to be curvature continuous are that there exist continuously differentiable functions \j and 
Kij and continuous functions [Hj and Vij, i = 1, j = 2, 4 and i = 2,4,j = 3, satisfying the equations 
A3.13i - D.16i and that the following relations 

2A4 j 3A' 1j4 K4 j 3 — Vx t 2 + ^4,3Al,4 + Ml,4«4,3 = 
2A 2) 3A 12 K 2j 3 — ^1,4 + ^2,3^1,2 + Ml.2^2,3 = 
2A 4j3 «;4 i 3K / 1i 4 - m )2 + ^4,3 + ^4,3«1,4 + ^1,4«4,3 = 
2A2,3K 2 ,3«l,2 - Ml,4 + ^2,3 + ^2,3^1,2 + ^1,2«|,3 = 
A 4)3 — A 2) 3A 12 + A4 j 3K / 14 — Al )2 K 23 + K 1A K 4,3 ~ Ml, 2 K 2,3 + ^1,4^4,3 = 
•^2,3 ~~ ^4,3Ai j4 + A2,3«l,2 ~~ ^1,4^4,3 + K l, 2^2,3 ~~ ^1,4^4,3 + ^1,2^2,3 = 

hold at the vertex V. 

Remark. We use here the same short notations as in Theorem 1, i.e., Ai )2 , ki j2 , Mi,2> ^1,2 etc, are to 
be interpreted as Ai^l), ki,2(1)> Mi, 2(1), ^1,2 (1)> etc > ie > as tne value of the functions at the vertex 
V. 

Proof The idea in the proof is to study the equation system consisting of (I3.11I >- (I3.16I > and to reduce, 
as far as possible, the number of equations including derivatives of second order. To simplify the 
notations in the proof we use r^ 1 ^ for r£ (u, v)\y = (1, 1), etc. 

First, we know that the equations (13.111 ) and (13.121 ) must hold. Using the relations (12.41 ) and (12.51 ) 



(3.17) 



in order to replace the vectors r$ and r^f 1 with r« and r„ respectively. Combining this with (12.21 ) 
and (12.31 ) in order to eliminate r^ 2 ^ and rffl we get 

A lj2 r« - r$ + A' li2 r« + + Ki )2 r& = 

A4, 3 rit ) - r[ 3 J + A 4i3 r« + <3(A li4 r« + rc lj4 r«) + Ki)3 r<$ = 

A lj4 r« - ri 4 j + 4 j4 r« + A' M r« + Kl , 4 r« = 

A 2 , 3 r$ - + 4 )3 (Ai )2 r« + /*i, 2 r«) + A' 2i3 r« + K 2i3 r$ = 0. 

In the above equation system we multiply the first equation by A 2; 3, the second one with —1, the third 
one with — A4 t 3 and add them to the fourth equation. We get 

(~^4,3 + ^2,3Ai j2 - A 4i 3k' m + Ai i2 K 2 3 - ki^k^) 1 "! + (A 23 + A 2i 3k' 1j2 - A 4i3 A / lj4 
+«1,2«4,3 ~ ^A K 4,3) r v^ ~ ^4,3«i,4rW + A 2j3 Ki i2 rW + K 2 , 3 r^ - «4,3r$ = 0. 
Considering the last part in the equation (13.181 ) and using the equations (13.131 ) and (13.161 ) it follows 
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that 

- X^Ki^r^ + A 2 ,3ti, 2 rl 1 t , ) + ^2,3ri 2 2 - K^r^ 
= -A4,3«i,4r^ + A 2 ,3«i,2r^ 

+K2,3(A? i2 r« + 2A 1)2 K 1 , 2 r[ 1 t ) + ^_ 2 r« + W , 2 r« + ^, 2 rW) 

-«4,3(Ai )4 rW + 2Ai )4 Ki,4r[ 1 t , ) + kI,4?uu + ^l^v + v i£ v u^) 
= «2,3(^i,2rl 1) + vx t 2T$) - K^m^r^p + ^1,4^^) 
+ ( K 2,3A? j2 - «4,3«1,4 - ^4,3«l,4)rW 
+ («2,3«1,2 - K4,3A? i4 + A 2j3 Ki )2 )r(, 1 t , ) 
+2(« 2 ,3Ai 2 «1,2 - «4,3Al,4«1.4)r^. 



In order to further reduce the above formula, we use the relations (12.61 ) and (12.71 ) to get 

K 2,3^1,2 — K4,3 K 1,4 ~~ A 4j 3Ki j4 = Kl,4Al )2 — «4,3«1,4 ~~ A4 5 3Ki j4 

= «l,4(Al,2 — A 4j 3 — Kl,4K 4 ,3) = 0, 

K 2,3«l,2 — K 4j 3Ai j4 + A2,3«l,2 = «2,3«1,2 ~~ K l,2^1,4 + A 2j 3Ki i2 

= «1,2(A2,3 - Al,4 + «1,2«2,3) = 

and 

«2,3Al,2«l,2 — «4,3Al,4«l,4 = «1,4«1,2 ~~ «1,2«1,4 = 0. 

Thus, it follows that 

- A4,3«i,4rS) + A 2) 3K ll2 r$ + ^2,3^ - «4, 3 rj$ 

= («2,3Ml,2 - K^l^)^ + (^2,3^1,2 - K^m)^. 

Input the above equality into the equation (13 - 1 8b - As before, the independence of the two vectors r„ 
and rjp gives 

— ^4,3 ~~ ^4,3^1,4 + A2,3Ai > 2 + Ai,2«2,3 ~~ K l,4«4,3 + ^2,3^1,2 _ ^4,3^1,4 = ^ 
^2,3 ~~ A 4j 3A / 1j4 + A 2 ,3^ / l,2 + w l,2«2,3 — ^1,4^4,3 + «2,3^1,2 ~ ^4,3/^1,4 = 0. 

The equations in (13.191 ) are necessary in order to fulfill the condition of geometric continuity of order 
2. 

Let us next study the equations (13- 1 3b — (13. 16b more closely. We rewrite two of these equations. 
Start with the equation (13.141 ) by adding to it the first equation in (13.111 ) multiplied by 2A4,3 M,3 and 
equation (13.131 ) multiplied by — A 4j 3/Ai j2 . We get the new equation 

(^4,3 + 2A4 ) 3/C4 )3 K / 1 4 - Ml,2T^)r^ + {2X^,3 ^,3^1 4 - u 1>2 t^)^ + ^4,3r[ 3) 

1,2 1,2 (3.20) 

+ (Ai 3 + 2A 4l 3«4,3«i,4 - A 1)2 A4,3)rW - 4 2 ^r« - (1 - ^)r® + = 0. 

M,2 M,2 

We continue with the equation (13.151) . To this one we add the first equation in (13.121) multiplied by 
2A 2j 3 k 2j 3 and equation (13.161 ) multiplied by — A 2; 3/Ai ;4 . We have 

(2^1,2^2,3^2,3 - ^1,4-^)^ + (/i2,3 + 2« / i i2 A 2 ,3«2 ) 3 ~ ^1,47^)^' + ^2,3^ 

M 1,4 (3.21) 

" «?,4^r« + (A 2 . 3 + 2 Kli2 A 2 ,3K 2 ,3 - A M A 2>3 )rW + ^r® - (1 - ^)r<? = 0. 
Ma Ma 
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Let us now consider the equation (13.20I ). We want to rewrite this equation in order to make it 
easier to handle. Using the relations (|231) - (I2.4I ). (I2.6M2.7I ) combined with (13.131 ) and (13.161 ). we get 

(^4,3^1,2 + 2Ai i2 A4 i 3K4 ) 3K / 1 4 - ^1,2X4,3)^ + A 4i 3(2Ai i2 K4 )3 Ai i 4 - Vl^T™ + ^4,3 Al, 2 r( 4) 

+ Ai i2 A 4 ,3(A4,3 + 2ki )4 K4,3 - Ai,2)rW - A4,3«i,2 r TO ~ (Ai,2 - A 4 ,3)rl 2 u ) + Ai^/^r^ 

= (^4,3^1,2 + 2Ai i2 A 4i 3K4 i 3K / 1) 4 - /Ui, 2 A 4j3 )r[ 1) + A 4j 3 (2Ai i2 K4 i3 A li4 - U lj2 )^ 

+ K4,3Ai,2(«i,4r£ L) + Ai^r^) + Ai )2 A4 ) 3« 1) 4«4 i 3r^ - A4 )3 rei )2 r$ - «i,4«4,3 (Ai )2 r^] 
+ 2Ai i2 Ki i2 rW + /e? )2 r$ + /x 1)2 rW + ^r^) + Ai^/^A^rW + 2A M rci, 4 rW 

+ <4,4^ul + MM 1 "^ + ^,4^) 
= (A*4,3Al,2 + 2Ai i2 A4 i3 K4 i 3k' 1j4 - / Ui i2 A 4i3 + 1^4,3 Al, 2 Kl,4 ~ «l,4«4,3Ml,2 + Ai^K^^m) 1 "^ 
+ (2Ai >2 A4 > 3K4 i 3A' 1j4 — Z^i j2 A4 j 3 + Z^Al^Ai^ — Kl,4K4,3Z/i i2 + Ai^K^Mi^Jr^ 
+ Ai >2 Ki i4 fi;4 i 3(A4 i 3 - Ai,2 + «1, 4^4,3)^2 - K? )2 (A4,3 + «1,4«4,3 ~ Ai j2 )r$ 

- 2Ai,2Ki,4«4,3(«l,2 - Ai,4/C4 )3 )rW 

= (A*4,3Ai, 2 + 2Ai i2 A4 i 3K4 ) 3K / 1j4 — /il )2 A 4i 3 + 1^4,3 Al,2^1,4 ~ «l,4«4,3Ml,2 + Ai^/cf^ 1,4)^ 
+ (2Ai i2 A 4i 3K4 i3 A 1j4 - 1/1,2 A4 )3 + ^4,3Al, 2 Ai i4 - « 1) 4K4 )3 I/i ) 2 + Ai^K^i^r^ = 0. 

The independence of the tangential vectors rl 1 ^ and r„ combined with the first relation in (12.71) gives 

Al, 2 (>4,3 + 2A 4i 3«:4 i 3K / 1i 4 - fJ, 1>2 + ^4,3^1,4 + ^1,4^,3) =0 ^ ^ ^ 

Al, 2 (2A4, 3 K 4 , 3 A / 1 4 — ^i j2 + Z/4,3Al,4 + ^1,4^4,3) = 0- 



We end our examination by simplifying the equation (13-21 b - Similarly as above, we use (12.2b . (12- 5b . 
d2~6T>— d277T>. combined with (|3~T3T ) and (l3~T6l We have 

A2,3(2A' 1i2 Ai >4 k; 2i 3 - z/i,4)rl X) + (M2,3Ai,4 + 2«/ 1)2 Ai,4A 2 , 3 «2,3 - Mi,4A 2 ,3)4 1) + "2,3^1,4*^ 

- A2,3«i,4 r w] + Ai,4A 2 , 3 (A 2 , 3 + 2ki, 2 «; 2) 3 - Ai, 4 )rW + \i,4,i4,3 r uu - (Ai,4 - A 2)3 )ry 

= A 2 , 3 (2A' 1 )2 A 1) 4/C 2 ,3 - ^l,i) r u + (A*2,3Ai,4 + 2k[ >2 Ai >4 A 2) 3K 2i 3 - /il,4A 2) 3)4 

+ ^2,3Ai,4(Ai, 2 r( t 1) + Ki, 2 rW) - A 2j3 ref 4 rW + A 1) 4A 2 ,3«i ) 2«2,3rE ) + *i,4«l,3(^i,2 r uu 
+ 2Ai )2 Ki j2 rW + K? i3 r$ + AH^r^ + ^r^) - rci,2«2,3 (Ai )4 r W + 2Ai i4 Ki,4r$ 
+ «l,4 r S u ) +W,4 r ^ 1) +^i,4ri 1) ) 

= (2A 2j 3A / 12 Ai j 4K 2j 3 — ^l,4A 2 ,3 + I/2,3Al,4Al,2 + ^1^1,4^2,3 ~ l / l,4«l,2«2,3)r^ 1 '' + (a*2,3Ai,4 
+ 2k/ 1j2 Ai )4 A 2 ,3K2,3 - Ml,4A 2 ,3 + f2,3Al,4«l,2 + ^1,2^1,4, K 2,3 ~ A*l,4«l,2«2,3)r^ + «i,4(Al,4 

- A 2j3 - /Cl,2«2,3)rS + Ai i4 Ki i2 K 2i 3(A 2i 3 + /Cl ) 2«2,3 - Ai,4)r£$ + 2Ai i 4/?1 )2 K 2i 3(k 2)3 Ai )2 

- «l,4)r^ 

= (2A 2 ,3A 12 Ai,4K 2 ,3 — ^l,4A 2 ,3 + J/ 2 ,3Al,4Al j2 + /X1, 2 Ai j4 k| )3 — l/l,4«l,2«2,3)ri + (A*2,3Ai,4 
+ 2k / 12 Ai i 4A 2i 3K 2i 3 - /Ul,4A 2 , 3 + Z/ 2) 3Al,4«l,2 + ^1^1,4^23 - Afl,4«l,2«2,3)4 = 0. 

Using the same argument as before, i.e., the independence of the tangential vectors rl 1 ^ and r^, 
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together with the second relation in (12.7b . it implies 

Al,4(2A2 ) 3A / 12 K2,3 — ^1,4 + ^2,3^1,2 + ^1,2^2,3) =0 
Al,4(/"2,3 + 2k' 1j2 A2,3K2,3 ~ A*l,4 + ^2,3^1,2 + ^1,2^3) = 0. 

Combining the results in (13.191) . (13.221 ) and (13.231 ) with the fact that 7^ 0, we get the compatibil- 
ity conditions (13.171 ) for G 2 , which are necessary and sufficient for having a simultaneous satisfaction 
of the equations (I3.11I) - (I3.16I ). This ends the proof. □ 



4 An algorithm 

In the previous part of this article we have achieved general results about regularity of 4-patch surfaces. 
In this section we will apply those results in the G 1 -case to a particular patch parametrization such 
as Bezier representation. Our goal is to make the connection between two patches more flexible, 
which will make it easier to handle geometrical and other demands. The way to handle this will 
be done by letting the functions A and k in (12.11) be suitable polynomials. We start this section by 
studying the consequences for the relations between the control points from a 2-patch surface in such 
a case. This result will then be used in creating a smooth 4-patch surface from an uncomplete such 
surface consisting of three patches. Here the compatibility conditions in Theorem[T]are of fundamental 
importance. The next step is to solve the problem of filling a hole in a surface in a smooth way, more 
precisely, we consider an uncomplete regular 9-patch surface as in Figure 5. To create the interior 
patch we use the same technique as in creating a 4-patch surface. Finally, what we have done so far 
will be used in creating a fillet surface, see Figure 6, which together with its surrounding becomes a 
smooth surface. 



r 




Figure 4: Two patches connected along a common boundary 



We begin by considering the well-known G 1 -relation (12.11 ) between two patches, see Figure|4] Let 
the two patches be given by (u, v) i-> r(u, v) and (u, v) (->• r(it, v), where u, v E [0, 1] and with the 
common boundary curve u i-> f(u, 0) = r(u, 1) for u E [0, 1]. Then 

f v (u,0) = \(u)t v (u, 1) + /c(it)r u (tf,l), 0<tt<l. (4.24) 

We restrict ourselves in using Bezier representation of the patches as follows 

r(n,v) = Elo E?=o %' * 4 (8) («) Bf (v) 

and 

*(«,«) = E to £7=o %■ B® (u) flj m) (v) , 
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where the parameters u, v € [0, 1] with m = 4,5, 6. The functions B^ for i = 0, 1, . . . , n, are the 
Bernstein polynomials of degree n, i.e., Bl n \u) = (™) (1 - «)*-* v? for u £ [0, 1]. 
Let the functions A and k be polynomials satisfying 

X(u) = A (1 - u) 2 + 2a (1 - u)u + A x u 2 

(4 25) 

K (u) = K (1 - u) 3 + 3/3 (1) (1 - u) 2 u + 3/3 (2) (1 - u)u 2 + ki u 3 , 

where Ao, \\, k q , k\ a, (3^ and f3^ are constants and the parameter u € [0, 1]. From formula (14.241 ) 
we get 



^(q;,! - q 4 , ) Bf\u) = (A (l - uf + 2a(l - u)u + \ lU 2 ) £ 3 (q^ - q i)2 ) B< 3) (t 

i=0 i=0 

2 

+ («o(l - n) 3 + 3/3« (1 - u) 2 u + 3/3^ (1 - u)u 2 + k x u 3 ) 3 (q m , 3 - q*, 3 ) sf 1 (« 



j=0 



= (Ao (1 - ~)(1 - \) (cu.3 - q il2 ) + 2a (1 - 1)1 ( qi _i, 3 - q^-i 



i=0 



+ Ai^^(q^ 2 ,3-q,-2,2)) J Bf ) (n) 

+ 3^ («o (1 - - - i) (q i+1 , 3 - q i>3 ) + 3/3W (1 - -1)1 (q 4 , 3 - <k_i^ 
i=0 

+ 3/?^ (1 - 5)4— (qi-l,3 " qi-2,3) + «i 5^ 3- (qi-2,3 " Qi-3,3)) Bf\u). 

Using the above identity we see that the control points must satisfy the next relations. We have 

m(qo,i - qo,o) =3 (A (qo,3 - qo,2) + « (qi, 3 - qo >3 )) 
™(qi,i - qi,o) =3 (^1,3 - qi,2) + ^ (qo, 3 - qo,2) 

+ — (q2,3 - qi, 3 ) H — 5— (qi,3 - qo,3jj 

m(q 2 ,i - q 2 ,o) =3 (-j^p (q2,3 - q 2)2 ) + (qi,3 - qi,2) + (qo,3 - qo,2) 

+ ^ (qs,3 - q 2 , 3 ) + — ^- (q 2 ,3 - qi,s) + — j^- (qi,3 - qo,3jj 



m(q 3 ,i - q 3 ,o) =3 (q 3 , 3 - q 3 ,2) + ^ (q 2 ,3 - q2, 2 ) + ^ (qi, 



(4.26) 



3-qi, 2 ; 



+ -jQ- (03,3 - q2,3j + "^q- (q2,3 - qi,3) + ^ (qi,3 ~ 00,3)) 



/- - \ /2qi . . 3 Ai . . 

^(q4,i - q4,o) =3 (— (q 3 ,3 - q3, 2 ) + -5— (q 2 ,3 - q 2 , 2 ) 

h — — (q3, 3 - q 2 ,3j + -5— (q 2 ,3 - qi, 3 )j 

m(q 5 ,l - qs.o) =3 (Ai (q 3)3 - q 3j2 ) + m (q 3;3 - q 2)3 )) • 
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Since the two patches r and r coincide along their common boundary curve we also have the 
following identity 

^q l ,0 J Bf ) (n) = ^q l , 3J Bf ) (n) 
i=0 i=0 

5 . . . . • - 1 

= E (t 1 " l^ 1 ~ \) *.3 + 2 (! " 5)7 *-i,3 + I— q-2,3) S f (B) («), u G [0, 1], 
i=0 

which, by identification, implies that 



qo,o = 


qo,3 






3 2 




qi,o = 


7 qi,3 + 7 qo,3 
5 5 






3 6 


1 


q2,o = 


10^3 + 7^1,3 + 


To 




1 6 


3 


q3,o = 


Tjj q3,3 + 75- qs.3 + 


To 




2 3 




q4,o = 


^q3,3 + ^q2,3 




q5,o = 


q3,3- 





(4.27) 

qi,3 



Thus, from the G l -condition between the two patches r and r, where r is known, the control points 
q^j, i = 0, 1, ... = 0, 1, are forced to satisfy the relations (14.261 ) and (14.27b - Nevertheless, first 
we have to decide the value of the parameters Ao, Ai, kq, k\, a, f3^ and . In the later applications 
it will be obvious how to choose certain of these parameters. 

In order to define the patch r completely, the other control points q^j, i = 0, 1, . . . , 5, j = 
2, ... ,m, must also be defined but not from the above regularity condition. There must be any other 
way to define them. 

The result we have achieved so far in this section will be used next, where we consider 4-patch 
surfaces. 



4.1 A 4-patch surface 

In the introduction of this section we considered a 2-patch surface. There we achieved relations 
between the control points of the two patches in (14.261 ) and (14.271) . We will use this when creating a 
4-patch surface. Looking at Figure |2j we assume that we have a G 1 -surface consisting of the three 
Bezier patches n 4 ), and r^ 2 ). Our goal is to create the fourth patch r( 3 ) in such a way that we get 
a G 1 -regular 4-patch surface. 

Since the patch r^ 3 ) must have a common boundary curve with r^ 2 ) and respectively, the 
second equation in (12.21 ) and (12.31) respectively must be fulfilled. We have 

if )(u,0) = A 2 , 3 (u)if)(u, 1) + K2 , 3 (n)if)(n, 1) 

r^(0,v) = X^(v)r^(hv) + ^ 3 (v)4 4 Hl,v), 

where the parameters u, v G [0, 1]. We continue with the next two equations which we get from the 
the second equation in (12.41 ) and (12.51) respectively. Those equations are the boundary conditions, i.e., 



r ( 3 )(u,0) = r( 2 )(u,l) 
A 3 \0,v) =r( 4 )(M). 



(4.29) 
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To be more specific, we exemplify by letting the patches r^, r^ 2 ) and r^ 4 ) all be Bezier patches of 
bi-degree (3, 3), while the patch r^ 3 ) is of bi-degree (5,5), i.e., 



5 5 

,(3) r( 5 )^ r(5) 



3 <"> 



r (3) (v) = EE< B n^ 

i=0 j=0 

where u, v € [0, 1]. 

Since the the patches n 1 ', r( 2 ) and r^ 4 ) have the same bi-degree the functions Ai i2 , Ai 4, rei 2 and 
Kl t 4, defined in the first equation of (12.21 ) and (12.31 ) respectively, must satisfy that Ai 2 and Ai t 4 are 
identically constant and and k± : 4 are polynomials of first degree. In fact these polynomials must 
satisfy Ki^(v) = «i 2(0)(1 — v) and Ki^ti) = ki^(0)(1 — u), but here we consider the case where 
the functions K12 and «i 4 are identically zero. There is in principle no difference between the two 
cases. 

Combining the above with the compatibility conditions (12.61 ) and (12.71 ) we get 



(4.30) 



•^2,3(0) — Ai 5 4 
^4,3(0) = Ai 5 2 
«2, 3 (0) =0 
«4, 3 (0) = 0. 

Let the functions A2, 3 , A4 ) 3, ^2,3 and K4 3 be defined as in (14.25b - Using (14.301 ) we have 

A 2 ,3(«) = Ai, 4 (1 - u) 2 + 2 a 2 ,3 (1 - u)u + A 2 , 3 (l) u 2 
A 4 ,3(v) = Ai, 2 (1 - v) 2 + 2 a 4j3 (1 - v)v + A 4)3 (l) v 2 
K2,z{u) = 30$ (1 - u) 2 u + 3$l (1 - u)u 2 + « 2)3 (1) u 3 

«4,3(«) = 3/3^ (1 - V) 2 V + 3$l (1 - V)V 2 + K4, 3 (l) « 3 , 



(4.31) 



where Aj^l), Ki, 3 (l) and for i = 2,4, are any constants, while the coefficients Oj 3 and 3 for 
% = 2, 4 need to fulfill certain compatibility conditions given in the end of this subsection, see (14.36b - 
The solutions of the system of differential equations (14.281 ) with boundary conditions (14.291 ) to- 
gether with the compatibility conditions (14.301 ) include the set of all possible patches that will result in 
a G 1 -regular 4-patch surface. In order to get such a solution r^ 3 ) we start by solving the first equation 
in (|4"29l From (I4T27T ) it follows that 



(3) 

q ,o = 


(2) 






Q (3) _ 
^1,0 — 


3 (2) 2 (2) 

5 q 1|3 + 5 qoi 






Q (3) _ 


10 q2 > 3 + 10 qi > 3 + 


1 

To 


(2) 

q ,3 


Q (3) " 


1 (2) 6 (2) 

10 q3 > 3 + 10 q2 > 3 + 


3 

To 


(2) 

qii 


Q (3) " 


2 (2) 3 (2) 

5 q3, 3 + 5 q2,3 






Q (3) _ 
^5,0 _ 


(2) 

q 3 , 3 - 







(4.32) 
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Considering next the first equation in (14.281 ) combined with (14.301) . it follows from (14.261) in the 
previous subsection that 

1o,l ~~ 1o,o = (loj ~~ loj) 

qg - qg -§ (qg - qg) + ^ (qg - qg) + ^ (qg - qg» 
qg - qg =| (qg - qg) + Mr <«3 - + <«8 - «S> 



+ ^(qg-qg) + ^P<qg-qg)) 



<& - <& =| (qg - qg) + ^ «8 - qg) + ^ c««8 - «S> (433) 



^ (qg - qg) + ^ (qg - qg) + ^ (qg - qg 



qg - qg -| (2f* (qg - qg) + ^aW (q g _ q «) 
+ ^(qg-qg) + ^(qg-qg)) 

„(3) „(3) _ 3 r\ MUn' 2 ' r^-U* flHn' 2 ) n( 2 )^ 
15,1 ~ 15,0 -g l A 2,3(lj (q 3 ,3 ~ 13,2^ + K 2,3(lJ (.1 3 , 3 ~ %,3JJ- 

When considering the second equation in (14.28b . using symmetry in (14.331) we immediately get the 
following relations 



Q (3) 
11,0 


(3) 

" 1o,o 


= 3 g 1>2 (1^-lg) 










_(3) 
11,1 


(3) 

- 1o,i 




2a 4,3 , (4) 

5 (q3 <° 


-lg) + 


3/3^3 ( 4 ) 


ig)) 


_(3) 
11,2 


(3) 
" 10,2 


=f( 3 £' a dS-iS) + 


6a 4 ,3 , (4) 
10 ^ 


-lS) + 


A 4,3(l) , (4) 

io lq3 -° 


-ig) 



+ ^#(qg-qg) + ! #(qg-qg)) 



qg - qg =| (qg - qg) + ^ (qg - qg) + ^ (qg - qg ) <■**> 

+ ^ (qSS - qg) + ^# (qg - qg ) + ^ (< $ - «$8>) 

qg - qg =| (qg - qg) + ^ (qg - qg) 

?(2) 



+ ^(qg-qg) + ^(qg-qg)) 
qg - qg =! (M3(D (qg - qg) + ^(i) (qg - qg))- 
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Finally, from the second equation in ( 14.291 ), similarly as in (14.321 ) combined with symmetry, we 

get 



(3) 

<lo,o = 


Q (4) 






(3) 
%,1 = 


3 (4) . 2 ( 4 ) 

g q 3) i + g qs,o 






(3) 

q j = 


io q3 - 2 + 10 ^ + 


i 

To 


(4) 

q 3 ,o 


(3) 


io q3 - 3 + 10 q3 ' 2 + 


3 

To 


a (4) 

q3,i 


(3) 


5 q3,3 + 5 q3,2 






(3) 

q j = 


(4) 







(4.35) 



(3) 

A conclusion from the result above is that the control points q) r for i = 2, . . . , 5, j = 0, 1, and 
i = 0, 1, j = 2, . . . , 5, are uniquely defined. For i, j = 0, 1, on the other hand, the control points 
qj^ are defined in two different ways above. Concerning qg there is obviously no problem. The 
fact that we have used the compatibility conditions (14.301 ) implies equality in the two different ways 

(3) (3) (3) 

defining q\ q and q [ . The last control point q^ [ must be chosen in such a way that the relations 

(14.281 ) are fulfilled and consequently it is uniquely defined. Therefore, let us calculate q[ 3 ] in terms 

°f {qg }f j=2 * n tne two different ways that this control point has been defined in (14.331 ) and (14.341 ) 
respectively, in order to see what conditions are needed for fulfilling uniqueness. First, it follows 

Q 2i3 =25(qg - qg - qg + qg) = 9A M (qg - qg) + 6a 2 , 3 (qg - qg) 

+ 9/43 (qS " qS) " 15AM (qg - qg) 
=9A M (qg - qg + qg - qg - (qg - qg)) 

+ (6a 2 , 3 - 15A M ) (qg - qg) + 9/3« (qg - qg) 
=(9A 1)4 + 9/3$) A lj2 (qg - qg) + (6a 2 , 3 - 6A M ) (qg - qg) 

- 9Ai 54 Ai i2 (qg - qg). 

Second, we consider the same expression as above given in (14.34b . We have 

<?4,3 =25(qJ5 - qg - qg + qg) = 9A 1)2 (qg - qg) + 6a 4 ,3 (qg - qg) 

+ y P4,3 (q 3 ,i - q 3 ,o) - 15A i,2 (q 3 ,o - ^2,0) 

-QX (n^ n ( 4 )4. n ( 4 ) „( 4 ) f„( 4 ) n/ 4 ^ 

-9Ai, 2 (q 31 - q 30 + q 30 - q 20 - (q 21 - q 20 Jj 

+ (6a 4 ,3 - 15A lj2 ) (qg - qg) + 9/?g (qg - qg) 
=(9Ai, 2 + 9/3$) A M (qg - qg) + (6a 4 , 3 - 6A 1)2 ) (qg - qg) 

- 9Ai )2 Ai i4 (qg - qg). 
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Finally, using (12.91 ) we get 

Q 2 ,3 - Qa,3 =(9Ai, 2 /5g - 6(a 4 ,3 - Ai, 2 )) (qg - qg>) 

+ (6(a 2 , 3 - A lj4 ) - 9A M /?g) (qS " <$) 
+ 9Ai, 4 Ai, 2 (qg - qg) - 9A lj2 A 1>4 (qg - qg) 
- 9Ai, 4 Ai, 2 (qg - qg) + 9A li2 A 1)4 (qg - qg) 
= (9Ai, 2 ^3 ) - 6(«4,3 - A 1)2 )) (qg - qg) 

+ (6(a 2)3 - A lj4 ) - 9A M /3«) (qg - qg) = 0. 

The fact that the vectors qg — qg and qg — qg are linearly independent implies that the control 
point qjfl is uniquely defined if and only if the coefficients 0^3, 3 for i = 2,4, satisfy 

3Ai ,2/3/3 = 2a 4 3 — 2Ai 2 

:j (4-36) 
3Ai j4 /3^ 3 ; = 2a 2i3 - 2Ai )4 . 

With the above criteria fulfilled the control points qW for i = 0, 1, . . . , 5, j = 0, 1, and for i = 0,1, 

j = 0, 1, . . . , 5, are uniquely defined by (14.32b — (14.35b - Concerning the other control points q^- , 

i,j = 2,3,4,5, there is no criteria given here how to choose these. Use any suitable method to 

(2) 

determine those control points. The same is true for the parameters /3^ 3 with i = 2,4, while there 
may be certain restrictions concerning A^l), ^,3(1) for i = 2, 4. 

In the case, referred to on page 13, where Ki )2 (0) ^ or ki, 4 (0) 7^ 0, then the coefficient relations 
in (14.361 ) are replaced by 

3Ai >2 ^ = 2 «4,3 - 2Ai >2 - Ai j2 k m (0) 
3Ai, 4 /?g = 2 «2,3 - 2A M - A M k 1j2 (0). 

Our goal so far has been to create a 4 th patch in a non-complete 4-patch surface in such a way that 
the new surface is G 1 -regular in spite of the fact that the functions A 2j 3 and A 4j 3 are non-constant. We 
have a further goal with this construction, which will be seen in the next subsection, but for the just 
mentioned purpose we can in fact choose r( 3 ) of bi-degree (4,4). Let 

aw = 5 (Mo) + Mi)), 

for i = 2, 4. Then we get the polynomials 

A 2 ,3(^) = Ai )4 (l -u) + A 2 , 3 (l)« 
A 4 ,3(w) = Ai )2 (1 - v) + \±fl{l)v 

«2,3(«) = 2/3 2 , 3 (1 ~ u)lL + K 2) 3(l) U 2 
«4,3(«) = 2^(1 -U)U + «4,3(1)« 2 - 
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In this case necessary and sufficient conditions in (14.361) are replaced by the next conditions 



which combined with (14.301 ) imply uniqueness of the control point j. On the other hand, this later 
solution is not flexible enough in the sense that it can not be used in the next application. This follows 
from (14.371) in the context of the following subsection. 

In the next subsection we will use the result we just received in order to fill a hole in a surface. 

4.2 Filling a hole in a surface 

In this subsection we consider a surface with an interior hole, more specific an incomplete 9-patch 
surface. By that we mean a surface like the one in Figure 5 where the part denoted by r^ 5 ) is not 
included in the surface. We assume that the incomplete 9-patch surface is G 1 -regular. We also assume 
that each of the original 8 patches are represented by a Bezier polynomial of bi-degree (3, 3) and that 
along the boundary curve between any two of the patches holds that the function A is constant and 
the function k = 0. See (14.24b . The problem to consider here is to create an interior patch r^ 5 ) with 
as low bi-degree as possible keeping the G 1 -regularity of the complete surface. This is the situation 
Sarraga considered in the paper lfl6l . With his assumption he needed bi-degree (6,6) for A 5 \ In this 
subsection we construct a patch r( 5 ) of bi-degree (5,5). 

In order to create a Bezier patch r^ 5 ) of bi-degree (5, 5) we start by considering it as part of the 
4-patch surface consisting of rW, r^ 4 ), r^ 5 ) and r^ 2 ). From the assumption above we have, as in the 
previous subsection, that the functions Ai j2 , Ai )4 are constant and ki )2 = fti )4 = 0. We know that 
there exist compatibility conditions which are necessary and sufficient in order to get a unique solution 

(5) 

for the control points : for i = 0, 1, . . . , 5, j = 0, 1, and i = 0, 1, j = 0, 1, . . . , 5. Next we do 
the same with r^ 5 \ A 8 \ and r^ 6 ^. We start by defining the missing patch r^ 5 ) as part of the first 
4-patch surface. We have to solve 



2A li2 /3 2i3 = A 4)3 (l) - Ai, 2 = A 4)3 (l) - A 4 , 3 (0) 
2A M /3 4 ,3 = A 2 , 3 (l) - A M = A 2 , 3 (l) - A 2)3 (0), 



(4.37) 




X 2 , 5 (v)r^(l,v) + K2j5 (v)r^(l,v), 



(4.38) 



and 



r^(u,0) =r( 2 )(u,l) 
r( 5 )(0,t>) =r( 4 )(l,w), 



(4.39) 



with 



A 4j5 (ti) = A 4)5 (0) (1 - u) 2 + 2 a 4i5 (1 - u)u + A 4j5 (l) u 2 

A 2 , 5 W = A 2 , 5 (0) (1 - v ) 2 + 2 a 2 , 5 (1 - v)v + A 2>5 (1) v 2 

k 4>5 (u) = K4, 5 (0) (1 - u) 3 + 30$ (1 - u) 2 u + 3/3$ (1 - u)u 2 + /t 4j5 (l) n 3 

«2,5(«) = «2, 5 (0) (1 - vf + 3/3$ (1 - v) 2 v + 3/3$ (1 - v)v 2 + k 2 , 5 (1) 



(4.40) 



where the parameters u, v € [0, 1]. 
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From the previous subsection we know that in order to get a unique solution for the control points 

(5) 

q\j for i = 0, 1, . . . , 5, j = 0, 1, and i = 0, 1, j = 0, 1, . . . , 5, necessary and sufficient conditions 
are 

A 4 , 5 (0) = Ai, 2 , A 2 , 5 (0) = Ai, 4 fAA ^ 

(4.41) 

«4,5(0) = 0, K2, 5 (0) = 

and 

3A M /?g = 2(a 2j5 -A M ) 
3Ai, 2 ^ =2(a 4 ,5-Ai, 2 ). 




Figure 5 : Filling a hole in a surface 



As in the previous subsection, compare (14.321 ) and (14.35b . we next get the solution of the boundary 
control points for the patch r( 5 ) 



(5) 

q ,o 


(4) 
= 4o,3 






(5) 

4 ,o = 


Q (2) 
43,0 




_(5) 
41,0 


- 3 (4) , 2 (4) 

— 5 4l,3 ^ 5 4o,3 






(5) 

4 ,i = 


3 Q (2) , 2 (2) 
5 43,1 ' 5 43,0 




(5) 
42,0 


- J- + J. a ( 4 ) + 
_ 10 42,3 + 10 4l,3 + 


1 

10 


(4) 
4o,3 


(5) 
4 ,2 = 


10 43,2 ^ 10 43,1 T 


10 43,0 


Q (5) 
43,0 


- J_ a ( 4 ) + J. Q ( 4 ) + 
~~ 10 43,3 ^ 10 42,3 ^ 


3 

10 


a (4) 

4l,3 


(5) 
4o,3 = 


10 43,3 t 10 43,2 ^ 


10 43,1 


(5) 
44,0 


- 2 (4) , 3 (4) 
~~ 5 43,3 ^ 5 42,3 






(5) 
40,4 = 


2 Q (2) , 3 (2) 
5 43,3 ' 5 43,2 




(5) 
45,0 


(4) 
= 43,3 






(5) 
4 ,5 = 


Q (2) 
43,3- 





(4.43) 
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We get, compare (14.331 ) and (14.341) . part of its interior control points 

(5) (5) _3Al£ f (4) _ (4 ) 
40,1 40,0 — r v4o,3 4o,2^ 



and 





-qfo 


4( 3A , 1,2 foil 

5 v 5 v 1>d 


"4© 


+ 2 7(q?i 


(4), . 3/3^ 5 } (4) 

- qn 9) H — (q 3 ■ 

■MJ , Z / r v A 1 , 


- 41)) 




(5) 
~ 42 


_3 ,3Ai, 2 , (4) 

- 5 l 1Q ^,3 


-qg) 


+ 6 i 4 "(qS 


qSI) + A ^ 8 (qil 


q$) 






+ 6 ^ (a {4) 

+ 1Q (4 2 ,3 




+ ^ 4 ' 5 fa (4) 


-qS3» 




(5) 
43,1 


(5) 

" 4 3 ,o 


= 3 ( Al ' 2 (q? 3 
5 v 10 VM3 ' 3 


43l) + 6 7(42l- 


qg) + 3 ^ 8 (iS- 








+ 3 ^ fa {4) 

+ 1Q (4 3 ,3 




+ 6 f (qg 






(5) 

4 4 ,i 


(5) 

- 4 4 ,o 


_ 3 /2«4,5 , (4) 
" 5 ( 5 (4 3 ,3 


-qg: 


1 + 3A ; 8 (qS 




-<$)) 


(5) 
45,1 


-a (5) 

45,0 


= 3A 5 7 ' 8 









(5) _ (5) _3A M , ( 2 )_ (2)x 
4l,0 qo,0 ~~ g 143,0 42,0^ 



og - oS5 =| <4 2 1 - qg) + ^ (43 - oS) + ^ <iS - <$)) 

4 - 45 -I - oS) + ^ (4!! - oS> + ^ (45 - <$) 



+ 5 #(45-45>+ 2 #<45-45)) 



(4.44) 



(445) 



qS5 - 4 -I (45 - oS> + ^ (45 - 45) + ^ (45 - 45) 

(2) 

lS3 - lg =| (45 - ig) + ^ - qg) + ^ (qg - qg)) 

qS-qS=^(qg-qg). 

where we also have used the fact that A4 ) s(l) = Ay t 8 and ^4,5(1) = as well as A2,s(l) = A 3i 6 and 
k 2i5 (1) = 0, which all follow from (14-5 1 b - 

Let us continue with creating the next 4-patch surface containing the 4 patches r^ 5 \ r^ s \ A 9 ^ 
and r( 6 \ see Figure 5. With the natural change of indices we have here the same equations as in 
(14-3 8b — (14.39b together with (14.40I ). Thus, using the following compatibility conditions 



^6,5(1) — Ag,8, Ag i 5(l) — Xgfi 
K6,5(l)=0, « 8 , 5 (1)=0, 



(4.46) 
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and 

?(2) 



3 A 9i6 /3 6 5 - 2 (a 8) 5 - A 9j6 ) 
3 A 9>8 /3g 2 g = 2 (q 6j5 — Ag^) , 



combined with A6,s(0) = X^,2, «6,s(0) = and As,5(0) = Xr^, Kg^(0) = from (14.5 D , we easily 
achieve the next control points, 
similar way as before. We have 



(5) 

achieve the next control points, q 4 - • for i = 0,1,..., 5, j = 4, 5, and i = 4, 5, j = 0, 1, . . . , 5, in a 



and 



(5) 


(6) 

= 1o,o 






(5) 
15,0 = 


(8) 

1o,o 




_(5) 
11,5 


- 3 (6) , 2 (6) 
~~ 5 1l,0 + 5 10,0 






a (5) _ 
15,1 _ 


3 Q (8) , 2 (8) 
5 10,1 ^ 5 10,0 




(5) 
12,5 


- A a (6) + A a (6) + 

— 10 12,0 + 10 1l,0 + 


1 

10 


(6) 

1o,o 


(5) 
15,2 = 


3. a (8) + A a ( § ) + 
10 10,2 ^ 10 1q,i ^ 


A a ( § ) 
10 lo,o 


_(5) 
13,5 


- J_ a (6) ■ A a ( 6 ) + 

— 10 13,0 ^ 10 12,0 ^ 


3 

10 


a (6) 

1l,0 


(5) 
15,3 = 


J_ a (8) , A a (8) + 
10 10,3 ^ 10 10,2 ^ 


A a ( § ) 
10 1o,i 


(5) 
14,5 


- 2 (6) , 3 (6) 
~~ 5l3,0 ^ 5l2,0 






Q (5) _ 
15,4 — 


2 a (8),3 (8) 
5 10,3 ^ 5 10,2 




Q (5) 
15,5 


-a (6) 

_ 13,0 






15,5 _ 


(8) 
10,3' 





1o,4 1o,5 ~~ r ' (lo,0 1o,l) 

q<5 - -| (oS - q?!) + ^ <«$ - qffi) + ^ (iS - <® ) 
qg - 4? =| (^f (qg - qg) + ^ - <£».) + ^ «S - «$) 

+ ^#(qg-qS) + ! #(q < S-q$)) 



qg - qg -| (<$ - qg) + ^ (q g - q g) + ^ (q c»> _ 
+ 5 #(qg-qg) + 5 #(qg-qH)) 

(2) 

qg - qg =| (q$ - qg> + ^ (qg - qg) + !$t «$ - 
qg-qg^tqg-qg) 



(4.48) 



(4.49) 
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and 



(5) (5) _3A M , (8)_ (8K 
44,0 45,0 ~~ g 140,0 Vll,0i 

45 - 41 -| (45 - q!5) + 2s» (q g> _ q <5) + Ml (q « _ 
qS3 - 45 -| <45 - 45) + ^ (45 - 45) + iif (45 - 45) 

+ ^(45-45) + 2 #(45-45>) 
45 - 45 =| (45 - qS3) + ^ (45 - 45) + ^ (45 - 45) 

+ ^#(45-45) + ! #(45-qg» 



(4.50) 



(5) 


(5) 


=|( 2 7(45-45n 


3A 9i6 ( 8 ) 

5 (q °' 2 


-qS5) + 


(5) 
44,5 


-a (5) 


_3A 9j6 ( 8 ) ( 8 ) 
— r 1^0,3 Mi,3/- 







3 ^8,5 / (8) (8) 



When collecting all the above control point we may have lost uniqueness of the doubly defined 

(5) 

control points for i = 4, 5, j = 0, 1, and i = 0, 1, j = 4,5. Considering a 4-patch surface 
around each of the other two vertices of the patch A 5 \ i.e., the two 4-patch surfaces constituting of 
r ( 2 ) 5 r (5) > r ( 6 ), r( 3 ) and r^ 4 ), r^ 7 \ n 8 ), r( 5 ) respectively, we have necessary and sufficient conditions 
for uniqueness of the just mentioned control points. These conditions are 

^4,5(1) = ^7,8, ^8,5(0) = A7 ; 4 
«4,5(1) = °> K 8, 5 (0) = 

^6,5(0) = A3 : 2, A2,5(l) = A 3j 6 

«6,b(0) = 0, K2, 5 (l) = 



(4.51) 



and 

3(2) 



3 A7_4 /3 4 5 — 2 (a8,5 — 
3^7,8/3sl = 2(a 4 ,5 - A 7j8 ), 



(4.52) 



3A 3 , 6 /3$ = 2 (a 2 , 5 - A 3 , 6 ) 

3A 3 , 2 /3g = 2(q 6i5 -A 3 , 2 ). 

First we need to fulfil the compatibility conditions (14.42b . (14.47b . (14.521 ) and (14.53b . where we have 
four degree of freedom for the parameters 0^5, ^."5 and pj?) with i = 2,4, 6, 8. After decided the 
value of the parameters, we have partly defined an interior patch A 5 \ where the control points are 
uniquely defined except for the undefined control points • with i,j = 2,3. Those points must be 
chosen in some way. One way to chose the undefined interior control points is by defining a Coons' 



(4.53) 
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patch from the boundary curves combined with their derivatives, see ifToll . or use the next definition 

Q (5) _ (5) , (5) _ (5) 
^2,2 — 42,1 + Ml, 2 4l,l 

Q (5) _ (5) (5)_ (5) 

43.2 - 43,1 + 4 4 ,2 44,1 (/ . 5 .. 

a (5) _ Q (5) , Q (5) _ Q (5) 1 ' 

42.3 — 42,4 + 4l,3 4l,4 

a (5) _ _(5) , (5) _ (5) 

43.3 — 43,4 + 44,3 44,4- 

Thus, we have proved that it is possible to represent the interior patch r^ 5 ) as a Bezier patch 
of bi-degree (5,5) in such a way that the complete 9-patch surface is G 1 -continuous. On the other 
hand, the representation of r( 5 ) with non-constant 5 -functions given above is very dependent of 
non-trivial Kj ^-functions, because, if, on the contrary, the ki ^-functions are identically zero, then the 
Aj^-functions must be identically constant. This follows trivially from the compatibility conditions. 

Suppose that the K^-functions are identically zero, but the A^-functions are non-constant, then 
the patch r^ 5 ) must be of at least bi-degree (6,6). See also Sarraga liToll . In this case A4 5 and A2,5 are 
defined as 

A 4 , 5 (n) = Ai, 2 (1 - uf + 3 C 1 - u f u + 3 "I? (! " «)« 2 + A 7,8 v? 
A 2 ,5(«) = Ai, 4 (1 - vf + 3 (1 - v) 2 v + 3q§ (1 - u)u 2 + A 3 , 6 v 3 . 
The compatibility conditions will now be 

a 4,5 ~~ A l,2 

a 2,5 — 'dis- 
continuing around the other vertices in the patch r^ 5 ) we also get the following conditions 

(2) \ (2) > (1) \ 

a 6,5 - A 9,8, "4,5 - A 7)8 , a 6 5 - A 3j2 

(2) \ (1) \ (2) > 

a 8,5 — A 9,6; "8,5 — A 7,4, «2,5 — A 3,6- 

The above conditions must be complemented with (14.41 b . (14.461 ) and (14.5 lb . With these restrictions 
we get the unique control points of r^ 5 ^ except for the undefined points q>^, i,j = 2,3,4. Those 
interior points can be defined through a Coons' patch as above or in a similar way as in (14.54b . i.e., 

„(5) (5) (5) (5) 
42,2 _ 42,1 + 4l,2 4l,l 

„(5) _ (5) (5) _ (5) 
44,2 _ 44,1 t q 5]2 45,1 

Q (5) (5) (5) (5) 

42.4 — 42,5 T 4l,4 4l,5 

Q (6) _ _(B) , (5) _ (5) 
44,4 _ 44,5 + 45,4 45,5 

4g = + |(qg + qgl) - KqSS + qg) = + + <*g) " + <$) 

qffi = + + q$) - + «$) = + |(qg + <$) - + <$) 
qg = qg + ^(qg + <$) - + «$) = + + «$) - + <$) 
qg = qg + ^(qg + <$) - + <$) = <*S + K^l + <$) - + 45) 



and 
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and finally 

(5) 1/ (5) . (5) . (5) . (5)x 1/ (5) . (5) . (5) . (5) N 
(i) 

Observe, the parameters a\ 5 for i = 2, 4, 6, 8, and j = 1,2, are here completely decided from the 
compatibility conditions. 

A further observation is that if the patch r^ 5 ) is of bi-degree (4,4) the compatibility conditions, see 
(14.371) . in this case are not flexible enough to solve our problem. 

We will use the main result of this subsection in the next part, where we will study a way to create 
fillet surfaces. 




u 



Figure 6: Creating a fillet 



4.3 Creating a fillet surface 

Consider two surfaces of regularity G 1 , each consisting of the patches r ( 1+3n ) for n = 0, 1, . . . , N — 1 
and r( 3+3n ) for n = 0, 1, . . . , N — 1. We want to connect these two surfaces by creating a fillet surface 
in such a way that altogether there will be one complete G 1 -surface. We assume that all the patches so 
far are of bi-degree (3,3) and that the connection between any two patches satisfies that the A-function 
is constant and the K-function is identically zero. See Figure 6. 

Our first step is to create the patches r ( 2+6n ) for n = 0, 1, ... , \{N + l)/2] - 1, where [N/2] 
denotes the integer part of N/2. Let the patch r ( 2+6rt ) be defined as a Bezier patch of bi-degree (3, 3) 
connecting r ( 1+6n ) and r ( 3+6n ) in a G 1 -regular way. Furthermore, there are no non-zero ft-function 
towards the two neighboring patches. In the next step we create the patches in between, i.e., r ( 5+6n ) 
for n = 0, 1, ... , [iV/2] — 1, as we did in the previous subsection. All together, this completes the 
construction of the fillet surface. 



Acknowledgement. I like to thank Roger Andersson for valuable discussions. 
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